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ABSTRACT 

We present new XMM-Newton observations of two luminous and high accretion-rate radio-quiet active 
galactic nuclei (AGNs) at z ^ 2. Together with archival X-ray and rest-frame optical spectra of three sources 
with similar properties as well as 25 moderate-luminosity radio-quiet AGNs at z < 0.5, we investigate, for the 
first time, the dependence of the hard ( ^ 2 keV) X-ray power-law photon index on the broad Hj3 emission-line 
width and on the accretion rate across ^ 3 orders of magnitude in AGN luminosity. Provided the accretion 
rates of the five luminous sources can be estimated by extrapolating the well-known broad-line region size- 
luminosity relation to high luminosities, we find that the photon indices of these sources, while consistent with 
those expected from their accretion rates, are significantly higher than expected from the widths of their Hj3 
lines. We argue that, within the limits of our sample, the hard-X-ray photon index depends primarily on the 
accretion rate. 

Subject headings: galaxies: active - galaxies: nuclei - X-rays: galaxies - quasars: emission lines - quasars: 
individual (Q 1 346 -036, HE 22 1 7 -28 1 8) 



1. WHAT DETERMINES THE X-RAY POWER-LAW SPECTRUM OF 
ACTIVE galactic NUCLEI? 

It is widely accepted that the hard ( ^ 2 keV) X-ray emis- 
sion from active galactic nuclei (AGNs) is primarily produced 
via unsaturated inverse Compton scattering of UV-soft-X-ray 
photons from the accretion disk by a corona of hot, likely 
thermal, electrons (e.g., Haardt & Maraschi 1991; Zdziarski, 
Poutanen, & Johnson 2000; Kawaguchi, Shimura, & Mi- 
neshige 2001). The emitted X-ray photon spectrum in the 
«2-100keV energy range is best described by a power-law 
of the form Ne °^ E^^, where F, the photon index, takes a 
fairly constant value of r^2, and is predicted to be only weakly 
sensitive to large changes in the electron temperature and the 
optical depth in the corona (e.g., Haardt & Maraschi 1991). 
Throughout this work, we refer to F only as the photon index 
in the hard ( ^ 2 keV) X-ray band. 

Nandra & Pounds (1994) found a (F)=1.95±0.15 for a 
sample of bright Seyfert 1 galaxies, confirming the mean 
predicted F and its small dispersion. Brandt, Mathur, & 
Elvis (1997) have studied the X-ray spectra of a sam- 
ple of nearby AGNs, including low-moderate luminosity 
Narrow-Line Seyfert 1 (NLSl) galaxies, defined as type 1 
AGNs having Hj8 full width at half-maximum intensity 
(FWHM) < 2000 km s"' (Osterbrock & Pogge 1985). The 
addition of NLSls, which were previously overlooked due to 
observational biases, broadened the range of F values, and 
allowed Brandt et al. (1997) to find that F is anticorrelated 
with FWHM(Hj3). This extended the analogous relation be- 
tween FWHM(Hj3) and the effective X-ray spectral slope in 
the 0.1-2.0keV band to harder X-rays (e.g.. Boiler, Brandt, & 
Fink 1996; Laor et al. 1997). The pseudo-thermal soft-X-ray 
excess radiation may be responsible for cooling the corona, 
thus steepening the X-ray spectrum (e.g.. Pounds, Done, & 
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Osborne 1995), and the pronounced soft excesses of NLSls 
(e.g., Puchnarewicz et al. 1995) may readily explain their 
steep X-ray spectra and the F-FWH]Vl(Hj3) anticorrelation. 

The remarkable dependence of the spectral shape of the 
X-ray emission, originating ^ 30 Rs from the central engine, 
on the width of a broad-emission line region (BELR) line, 
emitted at fulO'^Rs, was interpreted as a fundamental de- 
pendence of F on the accretion rate, since, as explained be- 
low, FWHM(Hj3) is considered an accretion-rate indicator in 
type 1 AGNs (e.g., Boroson & Green 1992; Brandt & Boiler 
1998). A high accretion rate would increase the disk tem- 
perature hence producing more soft-X-ray radiation (the "soft 
excess") and, at the same time, increase the Compton cooling 
of the corona and steepen the hard-X-ray power law. 

Recent X-ray studies of nearby (z <, 0.5) type 1 radio-quiet 
quasars (RQQs) have consistently shown that F ~ 2.0 ±0.5, 
and some of those studies confirmed the Brandt et al. (1997) 
F-FWHM(Hj3) relationship (e.g., Leighly 1999; Reeves & 
Turner 2000; Porquet et al. 2004; Piconcelli et al. 2005, 
hereafter P05; Brocksopp et al. 2006). Photon indices of 
F ^ 2.0 are also observed for 0.5 ^ z ^ 6 RQQs (e.g.. Reeves 
& Turner 2000; Page et al. 2005; Shemmer et al. 2005), how- 
ever, the spread in their F values seems smaller than that for 
nearby sources (see Fig. 3 of Shemmer et al. 2005). On the 
other hand, Grupe et al. (2006) report (F) = 2.21 ±0.52 for 
a sample of z > 4 RQQs, and suggest that the relatively steep 
X-ray spectra of such sources may be attributed to high accre- 
tion rates, in analogy with NLSls. 

In this Letter we test the hypothesis that the accretion rate 
largely determines the hard-X-ray spectral slope in AGNs, 
and argue that previous studies have not been able to break the 
degeneracy between the dependence of F on FWHM(Hj3) and 
the accretion rate, since highly luminous AGNs were left out 
of the analyses (e.g., Porquet et al. 2004). This is performed 
by investigating X-ray and Hj3 spectral-region data for a well- 
defined sample of 30 moderate-high-luminosity RQQs, se- 
lected for having high-quality X-ray and Hj3 -region spec- 
troscopy. Five of the sources are luminous RQQs at z 2, al- 
lowing, for the first time in this context, expansion of the AGN 
parameter space by 3 orders of magnitude in luminosity and 
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TABLE 1 
XMM-Newton OBSERVATION LOG 



Quasar 


RA (J2000.0) 


DEC (J2000.0) 






Observation 
Start Date 


Net Exposure Time (ks) / Source Counts 
MOSl MOS2 PN 


Q 1346-036 
HE 2217-2818 


13 48 44.08 
22 20 06.77 


-03 53 25.0 
-28 03 23.9 


2.370 
2.414 


2.47 
1.28 


2005 Jul 05 
2005 Oct 31 


10.4 / 90 10.4 / 144 9.4 / 520 
25.1 /732 25.1 /790 21.2/2968 



^Systemic redshift measured from the optical emission lines and obtained from S04. 
*^Neutral Galactic absorption column density in units of 10-*^ cm~- obtained from Dickey I 

black-hole mass (Mbh)- Two of the z - 2 RQQs, Q 1346-036 
and HE 2217-2818, were selected for XMM-Newton obser- 
vations from the Shemmer et al. (2004; hereafter S04) sample 
of luminous and high accretion-rate quasars based on their 
expected high X-ray fluxes in that sample. In §|2lwe present 
the new XMM-Newton observations, their reduction, and the 
X-ray spectral fitting of these two sources; their X-ray and 
optical properties are presented in §|3] In §|4]we discuss our 
results, and in particular, the dependence of F on the accre- 
tion rate for AGNs. Throughout this work we consider only 
radio-quiet AGNs to avoid any contribution from jet-related 
emission to the X-ray spectra. Luminosity distances are com- 
puted using the standard cosmological model with parameters 
Q.A = 0.7, D.M ^ 0.3, and Ho = 70km s"' Mpc"'. 

2. XMM-NEWTON OBSERVATIONS AND DATA REDUCTION 

Tabled gives a log of the XMM-Newton (Jansen et al. 
2001) imaging spectroscopic observations for Q 1346—036 
and HE 2217— 2818; the data were processed using standard 
XMM-Newton Science Analysis System^ v6.5.0 tasks. The 
event files of the observation of Q 1346—036 were filtered 
to remove '^35ks of flaring-activity periods; the net expo- 
sure times in Table[2 reflect the filtered data (time filtering 
was not required for the HE 2217— 2818 observation). The 
X-ray spectra of the quasars were extracted from the images 
of all three European Photon Imaging Camera (EPIC) detec- 
tors using apertures with radii of 30". Background regions 
were at least as large as the source regions. The spectra 
of Q 1346-036 (HE 2217-2818) were grouped with a min- 
imum of 10 (50) counts per bin. Joint spectral fitting of the 
data from all three EPIC detectors for each quasar was per- 
formed with XSPEC vl 1.3.2 (Arnaud 1996). We employed 
Galactic-absorbed power-law models at observed-frame en- 
ergies >0.6keV, corresponding to >2keV in the rest-frame 
of each quasar, where the underlying power-law hard-X-ray 
spectrum is less prone to contamination due to any potential 
soft excess emission and ionized absorption. The best-fit F 
values and power-law normalizations from these fits are given 
in Table|2j and the data, their joint, best-fit spectra, and resid- 
uals appear in Fig. [J 

We searched for intrinsic absorption in each quasar by 
jointly fitting the spectra with a Galactic-absorbed power law 
model including an intrinsic (redshifted) neutral-absorption 
component with solar abundances in the same energy range 
as before. None of the quasars shows significant intrinsic ab- 
sorption; upper limits on intrinsic A^h values appear in Ta- 
ble|2] Fig. n includes a F-A^h confidence-contour plot from 
this fitting for each quasar By applying F-tests between the 
models including intrinsic absorption and those that exclude 
it, we found that neither data set requires an intrinsic absorp- 
tion component. We also note that our XMM-Newton spec- 
tra show no indication of Compton-reflection components or 

^ http://xmm.esac.esa.int/sas. 



Lockman(1990). 

Fe Ka lines; such features are expected to be relatively weak, 
below our detection threshold, in the luminous sources under 
study in this work (e.g., the Fe Ka equivalent width is ex- 
pected to be < 50 eV; Page et al. 2004a). 

Finally, we searched for rapid (^1 hr timescale in the rest 
frame) X-ray variations in the XMM-Newton data of the two 
quasars by applying Kolmogorov-Smirnov (KS) tests to the 
lists of photon arrival times from the event files, but no signif- 
icant variations were detected. 

3. X-RAY AND OPTICAL PROPERTIES 

Basic X-ray and optical properties of Q 1346— 036 and 
HE2217— 2818 are given in Table|2] The luminosity at a 
rest-frame wavelength of 5 100 A is given in column (6), and 
FWHM(Hj3) is given in column (7); these data were obtained 
from S04. Column (8) gives /Wbh, determined as 



lO^Mp 



= 4.35 



vLv(5100A) 



1044 ergs s- 



0.7 r 



FWHM(HJ3) 



T 2 



lO^k] 



ms 



(1) 



and based on the recent reverberation-mapping results of Pe- 
terson et al. (2004) and the Kaspi et al. (2005) BELR size- 
luminosity relation (see also Kaspi et al. 2000). We note that 
this relation relies on a sample of AGNs having luminosities 
up to RslO^'' ergs s^^, and extrapolating it to higher luminosi- 
ties is somewhat uncertain; a reverberation-mapping effort is 
underway to test the validity of such extrapolations (S. Kaspi 
et al., in preparation). However, as Mbh is expected to scale 
with luminosity, for a given Hj3 width, and since in this work 
we perform only nonparametric statistical ranking tests, our 
results are not significantly sensitive to the precise coefficient 
values in this relation. Using equation Q, the accretion-rate 
ratios (col. [9]), iboi/i'Edd (where Lboi is the bolometric lumi- 
nosity; hereafter L/Ledd), are given by 



L/LEdd = 0.15/(L) 



vLv(5100A) 



1044 



ergss 



-1 



0.3 



FWH]Vl(Hj8) 
lO^ki] 



.-1 



(2) 

where we have employed equation (21) of Marconi et al. 
(2004) to obtain f{L), which is the luminosity-dependent 
bolometric correction to vLv(5100A). The bolometric cor- 
rection is / 6 — 8 (/ ~ 5) for the z < 0.5 (z ~ 2) sources in 
this work. 

The optical-X-ray spectral slopes in column (10) are de- 
fined as aox= log(/2kev//2500A)/log(^2kev/v2500A)' where 



/2kev and 



2500 A 



are the flux densities at 2keV and 2500 A, 
respectively. The aox values were derived using the photon in- 
dices and fluxes in columns (2) and (3), respectively, and the 
optical luminosities in column (6), assuming a UV continuum 
of the form /y °^ V"°^ (Vanden Berk et al. 2001). The 
values for our sources are consistent with the expected val- 
ues, given their optical luminosities (e.g., Steffen et al. 2006). 
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TABLE 2 

Best-Fit X-Ray Spectral Parameters and Optical Properties 













log vLv(5100A) 


FWHM(Hj3) 


logMsH 






Quasar 


r 


/v(l keV)" 


(d.o.f) 




(ergss"') 


(kms"') 


(Mr,) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


Q 1346-036 


2.02 ±0.17 


1.6±0.2 


37 (52) 


< 1.2 


46.9 


5100 


10.0 


0.3 


-1.69 


HE 2217-2818 


1.97 ±0.06 


3.5±0.2 


50 (50) 


<0.6 


47.2 


5200 


10.3 


0.3 


-1.69 



Note. — The be.st-fit photon index, normalization, and x~ were obtained from a Galactic-absorbed power-law model. Errors represent 90% confidence limits, taking one parameter 
of interest (A^^ = 2.71). 

^Power-law normalization for the pn data, taken from joint fitting of all three EPIC detectors with the Galactic-absorbed power-law model, and is given as the flux density at an 
observed-frame energy of 1 keV with units of 10"^^ ergs cm"- s"' Hz"' . 

^Intrinsic column density in units of 10"" cm"". Upper limits were computed with the intrinsically-absorbed power-law model with Galactic absoiption, and represent 90% confidence 
limits for each value. 




Observed-Frame Energy (keV) 

Fig. 1. — Data, best-fit spectra, and residuals for the XMM-Newton ob- 
servations of Q 1346-036 (top) and HE 2217-2818 (bottom). Open circles, 
filled squares, and open squares represent the EPIC pn, MOSl, and MOS2 
data, respectively. Solid lines represent the best-fit model for each spectrum, 
and the thick line marks the best-fit model for the pn data. The data were fitted 
with a Galactic-absorbed power-law model in the observed-frame >0.6keV 
band, that was extrapolated to 0.3 keV in the observed frame. The % residuals 
are in units of a with error bars of size 1. The inset in each panel shows 68%, 
90%, and 99% confidence contours for F and N\i, when the data are fitted 
with an additional neutral intrinsic-absorption component. 

4. IS F AN ACCRETION-RATE INDICATOR? 

In Fig.|2lwe plot T vs. vLv(5100A), FWHM(Hj3), Mbh, 
and L/Leaa for Q 1346-036 and HE 2217-2818. Fig.Elalso 
includes 28 unabsorbed Palomar Green (PG) RQQs (Schmidt 
& Green 1983) which have high-quality XMM-Newton and 
optical data; 25 of the quasars are at z < 0.5 and three, 
PG 1247+267, PG 1630+377, and PG 1634+706, are at 
1.5 — 2. Photon indices in the 2-12 keV rest-frame band 
for 27 of the PG quasars were obtained from Table 3 of P05; 
the photon index of PG 1247+267 was obtained from Page et 
al. (2004b). Optical data for the PG quasars were obtained 
from Neugebauer et al. (1987), Boroson & Green (1992), 
Nishihara et al. (1997), and Mcintosh et al. (1999), and pro- 
cessed with equations Q and Q. The fact that the X-ray and 
optical data are not contemporaneous may provide a source 
for scatter in Fig.|2] However, since most of the PG quasars 
do not exhibit significant variations in Y (e.g., George et al. 
2000) and Hj3 width (Kaspi et al. 2000), and their optical- 
continuum flux variations are typically at a level of ^ 50% 
over timescales of several years (Kaspi et al. 2000), the data 
in Fig.|2lare not expected to be strongly affected by variability 
(see also equations Q and O). 

We computed Spearman rank-correlation coefficients be- 
tween r and vLv(5100A), FWHM(Hj3), Mbh, and L/Ledd 



for the 25 low-luminosity, z < 0.5 sources. We found that, 
with the exception of vLv(5100A), the investigated prop- 
erties are significantly correlated with F (with 99.9% con- 
fidence), in agreement with Porquet et al. (2004) and P05. 
The F-Mbh and F-L/Lem correlations are largely a conse- 
quence of the r-FWHM(Hj3) correlation (Brandt et al. 1997), 
since both Mbh and L/Lem depend strongly on FWHM(Hj3), 
and the sample spans a relatively narrow luminosity range 
[vLv(5100A) -10'*'*-10'*'^ergss-i; see eqs. Q and Q, and 
Fig.|2f!]. The absence of a F-L correlation has been observed 
for larger AGN samples spanning broader luminosity ranges 
(e.g., Shemmer et al. 2005; VignaH et al. 2005). 

We repeated the above correlations, adding the five lumi- 
nous quasars at z 2 to the analysis. We found that F re- 
mains uncorrected with vLy (5 100 A), and that the F-Mbh 
correlation disappeared. This is a consequence of extend- 
ing the luminosity and Mbh ranges by ~ 3 orders of magni- 
tude, while our measured F values are typical of high-redshift 
quasars (e.g., Shemmer et al. 2005). The F-FWHM(Hj3) cor- 
relation remained significant with 99.9% confidence, but the 
(anti-)correlation coefficient dropped from 0.61 to 0.44. The 
F-L/Ledd correlation remained significant, with 99.9% confi- 
dence, and the correlation coefficient, 0.60, has not changed. 
This supports the hypothesis that the X-ray photon index de- 
pends primarily on L/Leaa- 

To test the hypothesis that F depends primarily 
on L/Leaa, we considered AGNs from Fig.|2j7 with 
3500<FWHM(Hj3)<6000km s^^ [which is the FWHM(H)3) 
interval of the five luminous z ^ 2 sources as well as that of 
the majority of the S04 quasars], and checked the significance 
of the deviation between the F values of the five luminous 
sources and the F values of the eight z < 0.5 PG quasars 
in that FWHM(Hj3) interval. A Mann- Whitney (MW) non- 
parametric rank test shows that the F distributions of the five 
high-luminosity quasars and the eight moderate-luminosity 
quasars are significantly different with 99.9% confidence. 
We also performed a KS test on the two F distributions and 
found that the F values of the two groups of quasars cannot 
be drawn from a single distribution with 99% confidence. 
On the other hand, MW and KS tests showed that, in the 
0.2 ;S L/Leaa ^ 0.5 range (representing the range of L/Lem 
values of the five luminous quasars as well as that of the 
majority of the S04 sources), the F values of the two groups 
of quasars are not significantly different (as can be clearly 
seen from Fig.|2j/). 

Compton-reflection contributions to the X-ray spectra of 
all 30 quasars are expected to be negligible given their 
moderate-high luminosities (e.g.. Page et al. 2004a, but 
see also Jimenez-Bailon et al. 2005). However, Page et 
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Fig. 2. — The hard-X-ray photon index vs. (a) vLv(5100A), (b) FWHM(H/5), (c) Mbh, and (d) L/Lem. Open circles mark PG quasars at z < 0.5 from the 
P05 sample; the NLSls in that sample are marked with filled circles. Diamonds mark luminous PG quasars at z ~ 1.5 — 2 with ai'chival XMM-Newtoji data. 
Q 1346-036 and HE 2217-2818 are marked with squares. 



al. (2004b) detected Compton reflection in PG 1247+267. 
We reanalyzed the X-ray spectrum of this source, and 
found that the hard-X-ray excess has no significant effect 
on r, whether fitted with or without a Compton-reflection 
component (consistent with Page et al. 2004b). We also 
searched for systematic differences between the analyses 
of Q 1346-036 and HE2217-2818 and the P05 sources 
by reanalyzing the XMM-Newton spectra of PG 1630+377, 
PG 1634+706, and the eight z < 0.5 PG quasars with 
3500<FWHM(Hj3)<6000km s"'. We reproduced the P05 
photon indices in the rest-frame 2-12 keV band with no sys- 
tematic deviations from their values. Additionally, fitting 
the data of Q 1346-036 and HE 2217-2818 from rest-frame 
2-12 keV (observed-frame 0.6-3.5 keV) did not alter our re- 
sults significantly. We note that for five z < 0.5 PG quasars 
with 3500<FWHM(Hj3)<6000km s^^ the rvalues fromP05 
are consistent with those from Porquet et al. (2004), and are 
somewhat lower than those from Brocksopp et al. (2006), due 
perhaps to analyses differences. 

Our results suggest that F depends primarily on L/Ledd^ as 
pointed out by e.g., Brandt & Bofler (1998) and Laor (2000). 
The key difference between this work and previous studies 
(e.g., Porquet et al. 2004; Wang, Watarai, & Mineshige 2004; 
Bian 2005) is that we expanded the quasar parameter space 



by adding highly luminous sources to the correlations (see 
Fig-Eh)- Such luminous quasars are found only atz ~1^, and 
their Hj3 lines, required for Mbh and L/Ledd determinations 
(see e.g., S04; Baskin & Laor 2005), are shifted into the ob- 
served near-infrared band. X-ray and H/3 spectral-region data 
of additional high-luminosity quasars are crucial to break the 
degeneracy between the r-FWHM(Hj3) and F-L/Ledd corre- 
lations and conclude that the X-ray power-law photon index 
depends on the accretion rate. Ultimately, a large enough in- 
ventory of X-ray and Hj3 spectral-region data for luminous, 
high-redshift quasars will allow testing of the hypothesis that 
such sources are analogous to NLSls (e.g., Grupe et al. 2006). 
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